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Motivated by recent experiments on Al nanoparticles, we have studied the effects of fixed electron 
number and small size in nanoscale superconductors, by applying the canonical BCS theory for 
the attractive Hubbard model in two and three dimensions. A negative "gap" in particles with an 
odd number of electrons as observed in the experiments is obtained in our canonical scheme. For 
particles with an even number of electrons, the energy gap exhibits shell structure as a function of 
electron density or system size in the weak-coupling regime: the gap is particularly large for "magic 
numbers" of electrons for a given system size or of atoms for a fixed electron density. The grand 
canonical BCS method essentially misses this feature. Possible experimental methods for observing 
such shell effects are discussed. 
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As the technology for fabricating ultrasmall metallic 
grains steadily improves, the typical sample dimensions 
are approaching molecular dimensions . The present 
accessibility to samples ranging from nanoscale to bulk 
has renewed interest in various features of the solid state 
that may or may not survive the excursion to ultrasmall 
dimensions. In superconducting Al samples, for exam- 
ple, the ability to distinguish even and odd numbers of 
electrons through tunneling experiments has called 
into question the use of the grand canonical ensemble to 
describe the electron pairing in these ultrasmall samples 
within a model with equal level spacings and within 
the attractive Hubbard model In these works the 
BCS theory of pairing was formulated within the canon- 
ical ensemble, following the early treatment of nuclei , 
and some rigorous "quality control" was provided by ex- 
act studies [|j7ljl^Jlltl. 

Within the attractive Hubbard model we have found 
two prominent features that emerged from the canonical 
BCS treatment, both of which were verified by the exact 
solutions. The first is the existence of "negative gaps" 
for odd electron number grains. By this we simply mean 
that a tunneling bias less than the charging energy would 
be required to tunnel an electron onto a grain with an 
odd number of electrons. The second is the existence of 
what were termed "super-even" electron numbers, where 
the tunneling bias required to tunnel an electron onto a 
grain with certain even numbers of electrons would be un- 
usually high. In this letter we investigate these features 
for various bandstructures in two and three dimensions, 
as might apply to Al, and briefly discuss some possible 
experiments to observe in particular the "super-even" ef- 
fect. 

We have adopted the attractive Hubbard model, whose 
specifics are well known. The additional feature we in- 
clude here is the possibility of using both periodic bound- 
ary conditions (PBC) as well as open boundary condi- 
tions (OBC), which are more appropriate for small sys- 
tems. Either of these is accomplished through a unitary 
transformation to a basis that diagonalizes the kinetic 



energy term. The BCS variational calculation is then 
performed with a wave function containing pairs of time- 
reversed states (n, |) and (n, ].) |jl^. The even and odd 
wave functions with j/-pairs are given by Q 
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with Ne — 2v and Ne — 2;/+!, respectively. The contour 
integral is on any counterclockwise path that encloses 
the origin. For odd iVg, the blocked state m is chosen 
so that it gives the lowest energy for a given coupling 
strength. Details with PBC were given previously 0, 
and with OBC they will be given elsewhere. We calculate 
the ground state energy for three systems with electron 
number A^e, A'e -|- 1, and A"e -b 2, and evaluate the energy 
gap by the formula Aa?^ = {En^-i — ^.E^^ + En^+i)/2. 

In the grand canonical ensemble, the number of elec- 
trons is fixed only on average. Thus we must solve the 
gap equation. 
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along with the number equation, 
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for gap parameters {A„} and chemical potential /i. 
Here, Vnn,mm is the transformed interaction, £"„ = 
{in — + is the quasiparticle energy and e„ — 



^ is the single-particle energy mod- 



ified with the Hartree term. The gap is given by Aq = 
min {En) for a finite size system, that is, with quantized 
energy levels {e„}. 

In the weak to intermediate coupling regime, the en- 
ergy gap as a function of electron number (or electron 
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density rie) roughly reflects the single-particle density of 
states (DOS). Thus in a simple cubic (SC) lattice in three 
dimensions (3D), the gap in the bulk limit is a smooth 
function of rie that increases from zero at zero density 
p3[ to a maximum value at half filling. In Fig. |^(a) 
we show the single-electron DOS for a SC lattice of 
iV = 16 X 16 X 16 = 4096 sites with PBC (solid curve). 
The result has been smoothed by convolution with a 
normalized Gaussian, and it is very similar to the bulk 
density of states shown by the dashed curve. Although 
N = 4096 is fairly large, the energy gap behaves quite 
differently from what we expect for a bulk system as a 
function of electron density for weak coupling. Results 
for the grand canonical BCS gap, Ao/i, are illustrated 
in Fig. 1(b) for \U\/t = 2 and (c) for \U\/t = 1. For 
\U\/t — 2 the overall scale of the gap Aq as a function 
of Ue resembles g{e) shown in Fig. 0(a). However, it 
has many fine structures; discontinuities at small density 
and cusps at larger density. This non-smooth behaviour 
is a result of the discrete density of states, i.e., quan- 
tized energy levels {e„} and their degeneracy in a finite 
size system. Such quantum structures of Aq turn out to 
be prominent for weaker coupling strengths, as seen for 
\U\/t = 1 in Fig. |l|(c): in this case there are discontinu- 
ities in the gap for the entire range of density. The 
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FIG. 1. (a) Density of states smoothed by Gaussian con- 
volution with width 7 = 0.7t for 3D TV = 16^ = 4096 sites 
(solid curve). The dashed curve is the bulk density of states, 
(b) Energy gap obtained by the grand canonical BCS, Ao/t, 
as a function of electron density Ue for A'^ — 16'^ = 4096 sites 
and \U\/t = 2. (c) Same as (b) but for \U\/t = 1. 



discontinuities or cusps in the gap arise from finite level 
spacings, while their positions as a function of Ue and the 
magnitude of the gap are determined by the degeneracy 
of levels, as will be explained in detail shortly. 
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3D SC:N. 8x8x8 = 512 
|U|/t-1.5 GCBCS A„ 
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FIG. 2. Energy gap Ajv^/t from the canonical BCS 
(crosses) and Ao/t from the grand canonical BCS (solid 
curve) as a function of electron density rie for 3D 
N = 8^ = 512 sites and \U\/t = 1.5 (upper frame): the 
canonical results are shown for only some representative den- 
sities. The number of single-particle states (spin degeneracy 
not included) as a function of rie (lower frame): each of the 
discrete levels is plotted at the density that corresponds to the 
closed-shell configuration up to that level in the zero-coupling 
limit. The canonical gap A]v^/t for weak coupling has jumps 
at these densities as seen in the upper figure, while the height 
of these jumps is determined by the energy spacing to the 
next level. 

In Fig. 1^ (upper frame) the gap A^r^ obtained by the 
canonical BCS is shown with crosses, along with the 
grand canonical gap Aq (solid curve), as a function of 
electron density rig, for a nanoscale system in weak cou- 
pling. The most obvious new feature is the "negative 
gap" , for systems with an odd number of electrons. As 
was already mentioned, this result has already been ob- 
served in small Al grains |^. For the even numbered 
grains note that most of the results shown follow the 
discontinuous, step- function- like behaviour of Aq. How- 
ever, anomalously high values occur at densities where 
the grand canonical result has discontinuities. These 
anomalies follow from the analogue of shell effects for 
a finite lattice of electrons. In the lower part of Fig. ^ we 
plot the number of single-particle states as a function of 
He for this system. Each of the discrete levels is plotted 
at the density that corresponds to electron number for 
filling all the levels up to that particular level ("closed- 
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shell" configuration and the height is the degeneracy 
of the level without the spin factor of two. It is clear that 
the densities where the canonical gap has a jump are the 
ones that correspond to the closed-shell configurations. 

In a closed-shell configuration, the occupation of levels 
is mainly driven by the kinetic energy. The cost required 
to occupy higher energy states exceeds the gain due to 
the increased interaction. Moreover, a careful examina- 
tion of the energy gain due to pairing reveals two distinct 
sources, a Hartree-like term, and the explicit BCS pair- 
ing term. The latter is small compared to the former, so 
the loss in energy reduction due to less mixing of states 
is indeed quite small. The Hartree-like term continues to 
play a role, however, which is why the value of the gap 
in the closed shell configurations is approximately equal 
to half the level spacing (i.e. one would have expected 
an additional pairing energy). The same physics occurs 
within the grand canonical ensemble , though the dis- 
continuity is the best these equations can do to account 
for the closed shell configurations. 



a function of single-particle energy for the entire band. 
The one for SC (top frame) for negative energy is the 
same as shown in Fig. ||, except now it is plotted vs. en- 
ergy so that the level spacings are clearly visible; these 
in turn determine the height of the jumps in the canon- 
ical gap. In fact for such a small system there are only 
two distinct level spacings in SC. This is why in the gap 
shown in Fig. |^, there are only two anomalously high 
values for the gap. In FCC (middle frame) there is no 
particle-hole symmetry and the degeneracy is more con- 
centrated near the top of the band. Compared with SC, 
the jumps at closed shell configurations will be more en- 
hanced by the larger level spacings and more frequent 
for smaller density; in addition the gap for open shells 
will be smaller (on average) because of less degeneracies. 
The BCC (bottom frame) has particle-hole symmetry as 
in SC, but the degeneracy is concentrated around zero 
energy. As in FCC there will be jumps more frequently 
at smaller densities, but near half-filling the gap will be 
continuous with the Fermi level at zero energy. 
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FIG. 3. The number of states (spin degeneracy not in- 
cluded) as a function of single-particle energy e for 3D 
A'^ = 8^ = 512 sites for different lattice structures; SC (upper 
frame), FCC (middle frame), and BCC (bottom frame). The 
SC and BCC level structures have particle-hole symmetry, 
while FCC does not. 

Different level structures result in different shell struc- 
tures in the gap. In Fig. |3| we show the number of states 
for iV = 8^ = 512 sites with PBC for SC, FCC (face- 
centred cubic) and BCC (body-centred cubic) lattices, as 



2D N = 8 X 8 = 64 : PBC vs. OBC 
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FIG. 4. The canonical energy gap Ajv^/t with PBC 
(crosses) and OBC (squares) and the grand canonical Ao/t 
with PBC (solid curve) and OBC (dashed curve) as a function 
of electron density Ue for 2D TV = 8^ = 64 sites and \U\/t = 2. 
With OBC, due to lower symmetry and hence less degeneracy 
of single-particle levels, the canonical gap tends to have more 
jumps, but with less height. 

The level structure also depends on the boundary con- 
dition. In Fig. ^ we illustrate this for a small 2D system, 
where we compare PBC and OBC. In Fig. ^the canonical 
gap Aat^ is plotted (for all densities) with crosses (PBC) 
and squares (OBC), and the grand canonical gap Aq is 
shown with solid (PBC) and dashed (OBC) curves, for 
\U\/t = 2. The SC lattice in 2D with PBC has a large 
degeneracy (a singularity in the bulk DOS) at zero single- 
particle energy. With OBC a relatively high degeneracy 
remains at zero energy, but for nonzero energy there are 
more levels with less degeneracies, because translational 
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symmetry is absent. For iV = 64 most of the levels are 
doubly degenerate, while some have no degeneracy. This 
is why in Fig. ^ the canonical gap with OBC has jumps 
more frequently (often periodic between iVe multiples and 
non-multiples of four) with less height than with PBC at 
smaller density. We note again that the canonical gap 
for even iVe for open shells and that for odd Nf. look 
symmetric about the x-axis (for both PBC and OBC). 
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FIG. 5. The canonical gap Ajv,/t (circles) for 2D TV = Ns 
sites for quarter filling rie = 0.5 as a function of 1/N; for 
Ns = 8, 10, 12, 14, 16, 18, 20, 22, 24, 26 and for \U\/t = 10 
(upper frame) and \U\/t = 1 (lower frame). The grand canon- 
ical gap Ao/t for Ns = 1000 is shown on the ordinate with a 
filled circle. 

We can also see the shell effects for a fixed electron den- 
sity by varying the lattice size. In Fig. ^the canonical gap 
An, is shown as a function of 1/N {N = N^) for a 2D SC 
lattice with PBC for quarter filling = 0.5 (circles), for 
strong coupling (upper frame) and weak coupling (bot- 
tom frame). The grand canonical gap in the bulk limit 
is indicated with a solid circle on the ordinate. In strong 
coupling the gap hardly depends on the number of sites. 
In weak coupling the gap exhibits strong size dependence 
as can be seen for \U\/t = 1. For Ns = 10, 14, 18, 22 
and 26, quarter filling corresponds to a closed-shell con- 
figuration and this can be seen clearly for Ng = 14, 18 
and 22 as the big jumps, which reflect the level spac- 
ing in each case. In contrast the Fermi level is open 
for Ns = 8, 12, 16, 20, 24. Interestingly the gap for open 
shells also changes as a function of size in a non-smooth 



way jij]. By far, however, the transition to the bulk 
regime is dominated by the oscillations of the magnitude 
of the gap between open and closed shell configurations. 

In summary, we have examined the tunneling gap for 
three dimensional ultrasmall superconducting grains, as 
a function of electron density, coupling strength, and sys- 
tem size. In weak coupling, shell effects are particularly 
prominent, and should be observable in very clean grains 
at low temperature. An ideal experimental arrangement 
would allow one to vary the electron density over a wide 
range. In this way one could observe the large modu- 
lation of the gap and identify "magic numbers" of elec- 
trons corresponding to the electron densities with anoma- 
lously large gap. However, in practice we anticipate that 
through the use of a gate electrode, one can vary the 
electron density only by a small amount (though large 
enough to see even/odd effects [^). Hence one will have 
to rely on ion-implanting a distribution of grain sizes, 
and thus make use of Fig. || to correlate gaps of different 
magnitude with different grain sizes. A systematic search 
should yield grain sizes whose electron number lies near 
a "magic number" so that tunneling a handful of elec- 
trons (one by one) onto the sample controlled by a gate 
electrode will allow one to observe large changes in the 
gap, as illustrated in Fig. ||. 
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